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Abstract 

An axial next-nearest-neighbor XY model is studied as a model of 
chiral liquid crystals which exhibit many ferro-, ferri- and antiferroelectric 
tilted smectic phases. Depending on the values of interaction parameters, 
this model exhibits Ising symmetric (i.e., the tilt directions of directors 
are parallel or anti parallel) phases or XY symmetric phases. Phases with 
each type-of-symmetry show the character of devil's staircase, which has 
been observed in experiments. 

1 Introduction 

Several chiral liquid crystal (LC) compounds have found to exhibit intermedi- 
ate subphases, SmCpj^ (also called SmC*), SmCpj2 (also called AF), and other 
ferrielectric phases (FIh, FIl and FI) between ferroelectric SmC* and antifer- 
roelectric SmC^ phases [1,2]. All these phases are tilted smectic phases and 
they have approximately the same tilt angle from the layer normal. Each phase 
consists of a helical stack of unit-set-of-layers (unit cells) ; the helix pitch is typ- 
ically on the order of several /um [3]. Structures of SmC^, SmCpjj^ and SmCpjj 
among those are well-known [2,4-7]: the unit cells characterizing SmC*, SmC^, 
SmCpj]^, and SmCpjj consist of 1, 2, and 3, 4 layers, respectively. The early 
experiments [1,2,8] have identified that the c-directors (the projection of a di- 
rector onto layer-plane) of neighboring layers are parallel or antiparallel. This 
experimental result indicates that the system can be described with an Ising 
model. In fact, it has shown that an axial next nearest neighbor Ising (ANNNI) 
model, which has competing interactions, exhibits infinitely many subphases 
and can reproduce the transition sequence shown in experiments [9, 10]. How- 
ever, the assumption of the Ising symmetry on these LC systems is unnatural. 
In practice, the c-directors of layers are proved to deviate slightly from the 
parallel-antiparallel configuration by later high-resolution experiments [4-7] . In 



this context, various phcnomcnological theories have been earried out to clarify 
the transition phenomena of the intermediate phases [12-16]. In these models a 
continuous rotational freedom of tilt direction is taken into account. To a greater 
or less extent, these models arc generalization of ANNNI model and the models 
are called axial next-nearest-neighbor XY (ANNNXY) models. Although the 
ANNNXY model exhibits infinitely many subphases, the system does not un- 
dergo discrete transitions between these subphases but continuously changes. 
By introducing an external aligning field on ANNNXY model, the transitions 
between subphases are shown to be recovered [11,17]. However such external 
field is artificial and its origin is ambiguous. 

As a more realistic model applicable to the antiferroelectric materials, the 
ANNNXY model with a biquadratic interaction has been introduced and in- 
troductory remarks on the crossover from XY character to Ising character are 
given [18]. In the present paper, we investigate the ANNNXY model with the 
biquadratic interaction and disclose general properties of this model. 

2 ANNNXY Model with a Biquadratic Interac- 
tion 

Following ref. [18] , we introduce a Hamiltonian of the ANNNXY model with a 
biquadratic interaction as 

H = Y.Hu (1) 
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where the denotes the azimuthal angle of a molecule i in l-th layer. The 
parameter J(> 0) is the interaction parameter between molecules belonging 
to the same layer. The summation in the first term of eq.(2) is done over all 
neighboring pairs in the l-th layer. The parameters Ji and J2 are the interaction 
parameters between molecules in the nearest neighboring layers and next nearest 
neighboring layers, respectively. Let Ji be positive without loss of generality 
since the transformation Ji — Ji followed by {t/?;^*"*} Wi^^ +ttI} leaves the 
above Hamiltonian invariant. In the following, we set J2 to be negative. The 
last term of eq.(2) is the biquadratic interaction with an interaction parameter 
K{> 0). This biquadratic interaction term stabilizes the parallel and antiparallel 
configuration of c-directors. This Hamiltonian can be viewed as a truncated 
Fourier approximation of an exact Hamiltonian. 

Our model turns to the ordinary ANNNXY model for vanishing value of K. 
In the limit oi K —y 00, it turns to the ANNNI model. Thus for sufficiently 
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large K the ordered phase has parallel and antiparallcl c-dircctors; we call such 
phases Ising symmetric phases. We call the other phases XY symmetric phases. 

In the present paper, we ignore the helical structure of the chiral smectic 
phases, and assume these phases as simple periodic repetitions of unit cells. 
Thus a smectic phase is characterized by a set of order parameters in the unit 
cell, and may be labeled with its wavenumber q. We note that the transforma- 
tion {ipi — > + vr/} results \n q ^ 1/2 — q. Thus, e.g., (j = 1/3 (SmCpj]^) and 
g = 1/2 (SmC^) phases under negative Ji correspond, respectively, to g = 1/6 
and q = phases under positive J\ . 

At zero temperature, there are three stable phases [18]: 5 = phase, q = \/A 
phase, and a phase Q in which the wave number changes continuously in the 
range Q < q < 1/A. The phase diagram for temperature T = is shown in the 
inset of Fig. 1. In the phase q = 1/4, irrespective of the value K, the azimuthal 
angles in a unit cell can be written as </3i = </52 = and ipz = ipi = it, i.e., 
this phase has Ising symmetry. The stable distribution of azimuthal angle in 
phase Q is ^5; = 2TTql, and the stable wave number is determined by interaction 
parameters as cos2iTq = J\/A{J2 + K). 



3 Theoretical Methods and Results 
3.1 Mean Field Approximation 

In this section, we investigate the stability of phases at finite temperature within 
the mean field approximation. We choose the order parameters of Z-th layer as 

Cl = (cosifl 
Si = (smipl 

C; = ^C0s2(pf^^ 

Si = (^sm2^f' 

where the angular brackets indicate a thermal average over all molecules in a 
single layer. We can assume, without loss of generality, that in Ising symmetric 
phases the c-directors are parallel to a;-axis and thus si = Si = 0. The mean 
fields due to these order parameters are 

^1= l3 [JZCI + Jl(c(_l + C(+l) + J2{ci-2 + Q+2)] , 

r]i= f3 [Jzsi + Ji{si-i + si+i) + J2{si-2 + S1+2)] , 

fii = f3K{Ci^, + Ci+,), ^' 
VI = PK{Si-i + Si+i), 

where we assumed that each molecule interacts z molecules in the same layer, 
one molecule in each of the two nearest neighboring layers, and one molecule in 
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each of the two next nearest neighboring layers. We can derive the mean field 
free energy per layer [18] using the above order parameters and mean fields as 



1 ^ 
r> ^ 



- {JZCI + Ji(q_i + Cj+i) + J2{ci-2 + Cl+2)}CI 

+ \ {Jzsi + Ji{si-i + + J2{si-2 + Sl+2)}SI 

+ \k{Ci.i + Ci+i)Ci + \k{Si.i + Si+i)Si 



- ^^T^T{^i,r]i, 111,1/1) 



(5) 



where p is a period (number of layers in a unit cell). The function r in the last 
term of eq.(5) is the one-molecule partition function defined as 

Ti^uVhl^U'^i) = / exp cos (yj + ?7( sin (yj + cos 2(yj + i/( sin 2(y?) (6) 
Jo 

Prom the derivatives of the free energy with respect to the mean fields, we obtain 
a set of self-consistent equations: 

Cl = -^lTlT{^l,r}l,IJ,l,Ul), 

si = In T{£,i,rii,iJ.i,iyi), 

(7) 

Cl = -^\nT{^i,r]i,i^i,ui), 

Si = ■^\nT{£,i,r]i,nun), 

where , rji, fii, and i/i in the right-hand sides are functions of order parameters 
defined in eqs.(4). One of the solutions of eqs.(7) which minimizes the free 
energy gives a set of order parameters of the thermodynamically stable state. 
The transition temperature between ordered and disordered phases is 

Tc = (l + Pi(Cc))T,^'', (8) 
where and Pi{z) are, respectively, 

P,{«) = h{z)/Hz), (10) 

with the n-th order modified Bessel functions /„ {() . The variable is a solution 
of 

C = 2i^Pi(C)/re, (11) 
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thus Tc is given by solving cqs.(8) and (11) sclf-consistcntly. The above critical 
temperature is derived [18], under the assumption that the transition is second- 
order, by expanding the free energy ^{ci, si,Ci, Si) around ci = si = 0, Ci = C^°^ 
and Si = S^^^K and letting the quadratic terms of the expansion vanish at 
critical temperature. The order-disorder transition temperature deviates from 
eq.(8) if the transition is first-order. Indeed, as we will see below, the transition 
temperature is slightly higher than Tc when K is comparable to T^. 

In the following, we restrict our consideration to a set of major ordered 
phases with wavenumbers q = 0, 1/8, 1/6, 1/5, 1/4, and a disordered phase. 
Since we disciiss the relative stability of these phases rather than infinitely many 
phases, the resulting phase diagrams should be viewed as rough approximations 
of correct phase diagrams. However, we can expect, from the analogy of ANNNI 
model, that other than the above six phases will occupy relatively small regions 
on the phase diagrams, and that we can safely ignore these minor phases. 

In order to reduce the number of order parameters, we make a reasonable 
assumption that the ordered phases are highly symmetric. In precise, when the 
period p is odd (i.e., p = 5 in this paper), we assume that a set of c-directors 
of a unit cell has mirror symmetry with respect to a plane perpendicular to the 
smectic layers; when p is even, we assume a two-fold rotational symmetry of the 
set of c-directors, in addition to the mirror symmetry. 

In the following calculation, we set J = Ji and z = We take | J2I as the unit 
of energy. Each ordered phase has a set of (reduced number of) self-consistent 
equations (7). We solved the sets of equations numerically, and determined the 
thermodynamic stable phase by comparing the free energy per molecule. 

3.2 Phase Diagrams 

At i^T = 0, our model is reduced to an ordinary ANNNXY model. Then the 
system transits through the phases of < g < 1/4 continuously [12-16] and thus 
the successive phase transitions are never observed. We show the phase diagram 
for JsT = in Fig. 1, in which we considered the relative stability of g = 0, 1/8, 
1/6, 1/5, and 1/4 phases. The phases with q = 1/8, 1/6, and 1/5 have XY 
symmetry, while q = and 1/4 phases have Ising symmetry. It seems there is 
a contradiction between K-Ji diagram (inset) and T-Ji diagram at ivT = and 
T = 0, i.e., stable regions of g = and 1/4 phases are overestimated in T-Ji 
diagram; such a contradiction, which can be observed in all the following phase 
diagrams, is due to our approximation in which infinitely many minor phases 
are ignored. 

Since the eff'ect of the biquadratic interaction stabilizes the Ising symmetric 
phases, g = and 1/4 phases dominate larger regions as K increases. The phase 
diagram for K = 0.2, Fig. 2, shows that the Ising symmetric phases spread from 
the lower-temperature region. The phases with g = l/8, g = l/6, and q = 1/5 
still have XY symmetry. 

When i^/1 J2I exceeds 0.5, the phases with < g < 1/4 vanish for any Ji at 
zero-temperature (see the inset of Fig. 1). The phases with g = 1/8, 1/6, and 
1/5 change qualitatively at this stage. In these phases, Ising symmetric phases 
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Figure 1: Phase diagram in the T-Ji plane for K = 0. Each fractional number 
indicates the wavenumber q of relatively stable phase, and D denotes a disor- 
dered phase. The shaded areas indicate the phases with XY symmetry, where 
non-vanishing si and Si are allowed; while the unshaded areas with q = and 
1/4 indicate the phases with Ising symmetry, where si = Si = 0. The inset is 
the phase diagram in the K-Ji plane for zero temperature. The wavenumber of 
the phase Q continuously changes from to 1/4. 




Figure 2: The same as in Fig. 1 but for K = 0.2. 
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Figure 3: The same as in Fig. 1 but for K = 1.0. The phases with q=l/Q and 
1/5 exhibit Ising symmetric phases. The broken line indicates the boimdary 
between XY symmetric phases and Ising symmetric phases. (Inset): The order 
parameter ratio s/cof q= 1/6 phase as a function of temperature at Ji/| J2I = 
1.8 around the XY- Ising transition temperature. 



appear at lower temperature region and the Ising symmetric regions expand as 
K increases. Figure 3 is a phase diagram for K = 1.0. In this phase diagram, we 
can clearly observe the cross-over from XY symmetric phases to Ising symmetric 
phases. 

Figure 4 shows the phase diagram for K = 3.0. The Ising symmetric phases 
are spread out almost the whole region. This figure also shows a difference from 
Figs. 1, 2 and 3 in the order-disorder transition temperature. The tempera- 
ture Tc (the order-disorder transition temperature under the assumption of the 
second-order transition) deviates from when K exceeds unity (see eqs.(8) 
and (11)); the deviation arises from the region where A;Brc/|"^2| < K. However 
the Tc does not give correct order-disorder transition temperature. In fact, the 
actual transition temperature is slightly higher than the calculated Tc in the 
region where K and feB^c/l J2I are comparable, as shown in Fig. 4. This fact 
indicates that the order-disorder transition is first-order, rather than second- 
order. Indeed, we observed that the order parameters changes discontinuously 
at the region where the order-disorder transition temperature is different from 
T 
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Figure 4: The same as in Fig. 3 but for K = 3.0. The XY symmetric phase 

remains only in q = 1/8 phase. The dotted line aroimd J1/IJ2I — 0.5 and 
kBT/\J2\ ~ = 3 is Tc defined in eq.(8). The order (g = l/4)-disorder transi- 
tion temperature (indicated by sohd fine) is sUghtly higher than the Tc and the 
transition is first-order in the region where K and fceTc/l J2I are comparable. 



4 Summary and Discussions 

Wc studied the transition behavior of ANNNXY model with biquadratic in- 
teraction which can be applicable to the successive phase transitions of anti- 
ferroelectric smectics; the biquadratic interaction is interpreted as the second 
Fourier component of a pair directional interaction. The crossover from the XY 
symmetry to the Ising one is identified for several representative phases. The 
Ising symmetric phase has been suggested to appear from the low temperature 
region as the strength of biquadratic interaction is increased [18]. At the XY 
symmetric phase with period p, the difference of tilt directions of directors in 
successive layers, ipi^i—ipi, is shown to be 2'jt/p, which indicates a simple helical 
structure. 

Our results show that the ANNNXY model with biquadratic interaction 
reproduce the phase sequence observed in the experiments. However, we may 
not expect the quantitative agreement with experiments. As an example, let us 
consider the LC molecules in g = 1/6 phase, which is essentially equivalent to 
g = 1/3 phase for negative Ji, and thus this phase corresponds to the 3-layer 
SmCpj]^. EUipsometric experiments have shown that c-directors of this phase 
do approximately but not exactly lie in a co-plain [6] . In precise, the azimuthal 
angles of molecules in a period are —6, S, and tt, with 5 < O.Srad (~ 30°). In 
terms of our model, this fact corresponds to the situation where |.s;/q| ^ 0.5. 
The inset of Fig. 3 shows a typical behavior of s; /c; in q = 1/6 phase at the 
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XY-Ising transition temperature. As shown in Fig. 3, the ratio si/ci is rather 
large(~ 1.7) at high temperature and approximately zero at low temperature; 
we cannot find any region where si/ci is comparable to the experimental results. 
For quantitative agreement of the present results with experimental ones, some 
additional modification to the present model is required. We may expect to 
obtain quantitatively reliable results by introducing a chiral interaction together 
with long range interactions [9,19,20], as done by Olson et al. [21] to explain 
the 5 in the framework of phenomenological free energy. 

The profile of phase diagram for Ising symmetric phases in Fig. 4 suggests 
apparently a character of devil's staircase [9, 10] and thus the system undergoes 
successive phase transitions. For phases of XY symmetry in Figs. 2 and 3, 
however, it is not clear from our present analysis whether the system exhibits 
the successive phase transitions or not. Analysis of the stability of the soliton 
excitation can be a powerful tool to certify the successive transitions between XY 
symmetric phases. In the ANNNXY model under the two-fold external field, the 
character of devil's staircase is certified even in the ground state [11]. However, 
in the present model, the phase changes continuously at the ground state for 
small biquadratic interaction as shown in the inset of Fig. 1. In this respect, the 
phase boundary of main phase should be studied at finite temperature, which 
will disclose the character of successive phase transitions. 

References 

[1] T. Isozaki, T. Fujikawa, H. Takezoe, A. Fukuda, T. Hagiwara, Y. Suzuki 
and I. Kawamura, Phys. Rev. B 48 (1993) 13439. 

[2] T. Matsumoto, A. Fukuda, M. Johno, Y. Motoyama, T. Yui, S. S. Seomun 
and M. Yamashita, J. Mater. Chem. 9 (1999) 2051. 

[3] L. S. Hirst, S. J. Watson, H. F. Gleeson, P. Cluzeau, P. Barois, R. Pindak, 
J. Pitn(^y, A. Cady, P. M. Johnson, C. C. Huang, A-M. Levelut, G. Srajcr, 
J. Polhnann, W. Caliebe, A. Seed, M. R. Herbert, J. W. Goodby, and M. 
Bird, Phys. Rev. E 65 (2002) 041705. 

[4] P. Mach, R. Pindak, A.-M. Levelut, P. Barois, H. T. Nguyen, C. C. Huang, 
and L. Furenlid, Phys. Rev. Lett. 81 (1998) 1015. 

[5] P. Mach, R. Pindak, A.-M. Levelut, P. Barois, H. T. Nguyen, H. Baltes, 
M. Bird, K. Toyne, A. Seed, J. W. Goodby, C. C. Huang, and L. Furenlid, 
Phys. Rev. E 60 (1999) 6793. 

[6] P. M. Johnson, D. A. Olson, S. Pankratz, T. Nguyen, J. Goodby, M. Hird, 
and C. C. Huang, Phys. Rev. Lett. 84 (2000) 4870. 

[7] A. Cady, J. A. Pitney, R. Pindak, L. S. Matkin, S. J. Watson, H. F. Gleeson, 
P. Cluzeau, P. Barois, A.-M. Levelut, W. CaHebe, J. W. Goodby, M. Hird, 
and C. C. Huang, Phys. Rev. E 64 (2001) 050702. 



9 



[8] A. Fukuda, Y. Takanishi, T. Isozaki, K. Ishikawa and H. Takezoe, J. Mater. 
Chem. 4 (1994) 997. 

[9] M. Yamashita and S. Miyazima, Ferroelectrics 148 (1993) 1. 

[10] M. Yamashita, Ferroelectrics 181 (1996) 201. 

[11] M. Yamashita and S. Takeno, J. Phys. Soc. Jpn. 68 (1999) 1493. 

[12] V. L. Lorman, A. A. Bulbitch and P. Toredano, Phys. Rev. E 49 (1994) 
1369. 

[13] M. Cepic and B. Zeks, Mol. Cryst. Liq. Cryst. 263 (1995) 61. 

[14] X. Y. Wang and P. L. Taylor, Phys. Rev. Lett. 76 (1996) 640. 

[15] A. Roy and N. V. Madhusudana, Europhys. Lett. 36 (1996) 22. 

[16] M. Skarabot, Mojca Ccpic, B. Zcks, R. Blinc, G. Heppke, A. V. Kityk, and 
L Musevic, Phys. Rev. E 58 (1998) 575. 

[17] S. Tanaka and M. Yamashita, Ferroelectrics 245 (2000) 209. 

[18] M. Yamashita and S. Tanaka, Ferroelectrics 245 (2000) 217. 

[19] M. Yamashita, J. Phys. Soc. Jpn. 65 (1996) 2122. 

[20] M. Cepic and B. Zeks, Phys. Rev. Lett. 87 (2001) 85501. 

[21] D. A. Olson, X. F. Han, A. Cady, and C. C. Huang, Phys. Rev. E 66 (2002) 
021702. 



10 



